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Abstract. This paper is part of the series presenting the final results obtained by the ESO Imaging Survey (EIS) project. It 
presents new / and K s data obtained from observations conducted at the ESO 3.5m New Technology Telescope (NTT) using 
the SOFI camera. These data were taken as part of the Deep Public Survey (DPS) carried out by the ESO Imaging Survey 
program, significantly extending the earlier optical/infrared EIS-DEEP survey presented in a previous paper of this series. The 
DPS-IR survey comprises two observing strategies: shallow K s observations providing nearly full coverage of pointings with 
complementary multi-band (in general UBVRI) optical data obtained using ESO's wide-field imager (WFI) and deeper J and K s 
observations of the central parts of these fields. Currently, the DPS-IR survey provides a coverage of roughly 2. 1 square degrees 
(~ 300 SOFI pointings) in K s with 0.63 square degrees to fainter magnitudes and also covered in J, over three independent 
regions of the sky. The goal of the present paper is to briefly describe the observations, the data reduction procedures, and to 
present the final survey products which include fully calibrated pixel-maps and catalogs extracted from them. The astrometric 
solution with an estimated accuracy of < 0"15 is based on the USNO catalog and limited only by the accuracy of the reference 
catalog. The final stacked images presented here number 89 and 272, in J and K s , respectively, the latter reflecting the larger 
surveyed area. The J and K s images were taken with a median seeing of 0"77 and 0!'8. The images reach a median 5o" limiting 
magnitude of J^g ~ 23.06 as measured within an aperture of 2", while the corresponding limiting magnitude in K s ^b is ~ 21 .41 
and ~22.16 mag for the shallow and deep strategies. Although some spatial variation due to varying observing conditions is 
observed, overall the observed limiting magnitudes are consistent with those originally proposed. The quality of the data has 
been assessed by comparing the measured magnitude of sources at the bright end directly with those reported by the 2MASS 
survey and at the faint end by comparing the counts of galaxies and stars with those of other surveys to comparable depth and to 
model predictions. The final science-grade catalogs together with the astrometrically and photometrically calibrated co-added 
images are available at CDS***. 
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1. Introduction 

Deep multi-wavelength observations of selected regions of the 
sky from space and ground-based facilities, combined with 
spectroscopic observations from large-aperture telescopes, of- 
fer the most promising means to probe the distant universe and 
study in a comprehensive way the evolution of galaxies and 
large-scale structures over a broad interval of look-back time. 
Foreseeing the need to complement observations being carried 
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out with the HST, Chandra, XMM, Spitzer and GALEX tele- 
scopes with ground-based multi-passband data in optical and 
infrared, the Working Group for public surveys (SWG) at ESO 
recommended the ESO Imaging Survey (EIS, Renzini & da 
Costa, 1997) project to undertake deep, optical/infrared obser- 
vations. 

This paper is part of a series describing recent data releases 
from the ESO Imaging Sur vey using the n ewly implemented 
EIS Data Reduction System (Ida Costa et alll2004l) . The reduc- 
tion of op tical data and details of the released products are de- 
scribed in iDietrich et all f 2006. hereafter Paper I) in the context 
of the optical foll ow-up of the XMM-Newton Serendipitous 
Sky Survey. In lOlsen et alJ J2006, hereafter Paper II) the reduc- 
tion of infrared data is described based on the infrared part of 
the EIS-DEEP survey. The last survey discussed in the present 
series is the EIS Deep Public Survey which combines optical 
( Mig nano et all E0O6I hereafter Paper III) and infrared obser- 
vations, with the latter being the topic of the present paper. 

The original infrared observations, conducted in 1998, tar- 
geted the HDF-S and CDF-S regions and fully calibrated im- 
ages, source catalogs and high-redshift galaxy candidates were 
immediately made public prior to the S cience Verification of 
the first unit of the VLT (see preprints by Ida Costa et allll998l 
iRengelink et all [l998). In spite of some shortcomings in the 
original reduction of the infrared images (see Paper II for a 
more detailed discussion) these were used in combination with 
the optical data, taken with the SUSI2 detector, to produce 
color catalogs and to color-select high-redshift galaxy candi- 
dates. These candidates were later used in spectroscopic obser- 
vations, conducted durin g commissioning of th e VLT and the 
FORS spectrograph (e.g., ICristiani et all 12000). which, in gen- 
eral, confirmed their estimated photometric redshifts. However, 
the total area of the survey was small and, in particular, only a 
small fraction of the field of view of Chandra was covered by 
the available data. 

In order to expand the original infrared coverage of EIS, 
complementing other ongoing optical surveys, and provide full 
coverage of the CDF-S < 350 square arcmin field and of the 
flanking regions, the SWG designed the following strategy for 
the Deep Public Survey (DPS). The EIS DPS survey consists of 
two parts: first, a deep, optical multi-passband (UBVRI) survey 
reaching limiting magnitudes iiiab ~ 25 mag and covering three 
distinct stripes (DEEP1, DEEP2 and DEEP3) of the sky; sec- 
ond, the contiguous coverage in infrared of the same regions, 
the focus of the present paper. The infrared coverage also has 
two parts: regions covered in K s down to a proposed 5o lim- 
iting magnitude of 21.3 (in the AB system) as measured in a 
2"aperture, and a smaller area covered down to Jab = 23.4 and 
K s ab = 22.7 mag. As explained in Paper III the three regions 
were selected for the following reasons: 1) DEEP1 because it 
overlaps with a deep radio survey; 2) DEEP2 because it in- 
cludes the CDF-S field (DEEP2c); and 3) DEEP3 due to its 
convenient characteristics and location in the northern galactic 
cap. The optical observations were carried out using the wide- 
field imager (WFI), mounted on the ESO/MPG 2.2m telescope 
at La Silla covering about 0.25 square degrees per pointing. 
Each stripe consists of four adjacent pointings (denoted by a-d 
in decreasing order of right ascension), yielding an area of one 



square degree. Further details about DPS can be found in the 
EIS web pages 1 . 

The primary goal of the DPS survey is to produce a data 
set from which statistical samples of galaxies can be drawn 
to study the large scale structures at high redshift. However, 
these data are also valuable for many other areas of inves- 
tigation, in particular for the cross-identification with X-ray 
sources detected by Chandra and XMM, as a photometric ref- 
erence for the ISAAC mosaic built by the GOODS s urvey 
llDickinson A GiavaW Eool iGiavalisco et all Eool. and 
combined with the optical data described in Paper III to im- 
prove photometric redshifts and to search for rare high-redshift 
quasars and galaxies (e.g., Kon getallhood) . 

The purpose of the present paper is to describe and present 
the data from new infrared observations carried out as part of 
DPS; the optical data are presented in a separate paper of this 
series (Paper III). The shallow K s survey is one of the largest 
currently available, and the d eep JK S is the largest to similar 
depth (see Elston et all 120061 for a review). The paper is or- 
ganized as follows: Sect, previews the overall observing strat- 
egy and describes the observations. Sect.|3]briefly describes the 
techniques employed to process the images, and to astrometri- 
cally and photometrically calibrate them. This section also in- 
cludes a discussion of the assessment made of the quality of the 
images and of the photometric calibration. Sect, ^presents the 
final products which include the stacked images and the single- 
passband catalogs extracted from them. Even though the obser- 
vations are planned as mosaics, the creation of image mosaics 
and associated catalogs are beyond the scope of the present pa- 
per. Also left for future papers are the production and analysis 
of optical/infrared color catalogs. While the goal of this paper 
is not to interpret the data, Sect, [^presents the results of com- 
parisons between the present observations and those of other 
authors. This is done to assess the quality of the astrometry and 
photometry of the data set. A brief summary is presented in 
Sect.il 

2. Observations 

The near-infra red part of the DPS wa s carried out using the 
SOFI camera ( Mo orwood et all 1998) mounted on the New 
Technology Telescope (NTT) at La Silla. SOFI is equipped 
with a Rockwell 1024 2 detector that, when used together with 
its large field objective, provides images with a pixel scale of 
0"29, and a field of view of about 4'.9 x 4'.9. The infrared point- 
ings were planned to produce a contiguous coverage of the re- 
gions of interest with sufficient overlap between adjacent point- 
ings to enable the construction of an image mosaic. The size 
of the overlaps was chosen to avoid a significant decrease in 
the effective exposure time due to the jitter pattern used, and 
to enable a common photometric zero-point to be established. 
As mentioned earlier the main objective of the survey was to 
cover each WFI pointing with shallow (720 sec) K s observa- 
tions, complemented by 3600 sec exposures in J and longer 
(4800 sec) exposures in K s over a 4 x 4 (~ 400 square arcmin) 
mosaic of the central region. 
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Field 


Ra (J2000) 


Dec (J2000) 


^-shallow 


J 


^s-deep 


Deep la 


22:55:00.0 


-40:13:00 


49 


16 


16 


Deep lb 


22:52:07.1 


-40:13:00 


49 


16 


16 


Deep2a 


03:37:27.5 


-27:48:46 


29 






Deep2b 


03:34:58.2 


-27:48:46 


49 


16 


16 


Deep2c 


03:32:29.0 


-27:48:46 


49 


12 


12 


Deep3a 


11:24:50.0 


-21:42:00 


36 


16 


16 


Deep3b 


11:22:27.9 


-21:42:00 


36 


16 


16 



The data presented in this paper were obtained as part of 
the ESO Large Programmes: 164.O-0561 and 169.A-725, cov- 
ering a total of 8 ESO periods (64-73). Data for this survey 
were taken in 12 visitor mode runs, comprising 81 nights in 
the period from January 31, 2000 to September 8, 2004. The 
observations include 14567 science exposures, out of a total of 
over 27000 frames, corresponding to a total of about 243 hours 
on-target. A full log of the observations can be found in the EIS 
web pages. 

The observations were carried out in groups of exposures, 
referred to as observation blocks (OBs), adopting the com- 
monly used jitter strategy as implemented in the ESO standard 
observing template (AutoJitter) for the SOFI instrument. Using 
this template, offsets are generated randomly within a square 
box of a specified side length, chosen originally to be 40", ap- 
proximately 15% of the SOFI detector field of view, which was 
later increased to 70"(starting November 2000) to improve the 
background estimate. A total of 733 science OBs, ranging from 
12 to 60 exposures, with individual exposures from 10 to 60 
seconds each, and 1169 photometric standard OBs typically 
consisting of 5 exposures (for a total of 100 sec) have been 
executed as part of this program. 

The field centers and final status of the observations are pre- 
sented in Tabled The table lists in Col. 1 field name, in Col. 2 
and 3 the position of the center of the field in J2000, and in 
Cols. 4-6 the number of pointings covered for each strategy. 
Typically the distance between adjacent pointings is 4'5 and 
the effective area of each pointing covers about 20-25 square 
arcmin. The numbers listed in the table represent the num- 
ber of observed pointings regardless of their exposure times. 
Completeness maps for each combination of strategy and field, 
showing the ratio between actual and planned integration times, 
can be found on the EIS web pages. Note that for the Deepl 
fields, a number of pointings reach only 85% of the planned in- 
tegration time for the K s deep-strategy. In the case of Deep2c, 
only 12 pointings are listed in the table because the remaining 
four are part of the EIS-DEEP survey (see Paper II). In keep- 
ing track of the layout of the mosaics, the case of Deep2c is 
somewhat unique as the larger mosaic was built from a smaller, 
pre-existing one. 

From the table one finds that there are 297 pointings 
contributing to the K s shallow strategy, corresponding to 
~ 2.1 square degrees, neglecting overlaps between pointings. 
In addition, there are 92 pointings in both / and K s contributing 
to the deep strategy. This corresponds to ~ 0.63 square degrees. 
Note that in the case of the K s observations the deep pointings 



have been included in the area estimate of both shallow and 
deep strategies. 

In Sect. [4] the final tally of the actual coverage achieved 
by the survey will be presented after taking into account the 
quality of the reduced data and the final images produced by 
the system. 

3. Data reduction 

The data reduction was carried out using the unsupervised 
mode of the EIS Data Reduction System ( Ida Costa et all 
2004), to reduce nightly images, to determine photometric so- 
lutions on a nightly basis, and to produce final stacked images 
and science-grade catalogs. A full description of the methods 
used in processing the images can be found in Paper I, and 
an extensive discussion on the validation of the software in 
handling infrared data is presented in Paper II. The interested 
reader is referred to these papers for more details. 

The astrometric calibration of the infrared frames was per- 
formed using the USN O-B catalog as refere nce following the 
method developed by Dia mdii et alJ dl993h . The method is 
based on a multi-resolution decomposition of images using 
wavelet transforms. This package has proven to be efficient and 
robust for pipeline reductions. The estimated internal accuracy 
of this technique is about half a pixel (~ 0"15) and is limited 
by the internal accuracy of the reference catalog. 

The photometric calibration for each night is based on ob - 
servations of standard stars taken from IPersson et alJ 11998). 
Typically two to seven stars were observed over a range of air- 
mass. For all nights independent photometric solutions were 
obtained using the photometric pipeline of the EIS data re- 
duction system. Depending on the airmass and color cover- 
age, linear fits with one to three parameters are attempted. If 
no standard star observations are available a default solution is 
assigned based on the solutions present in the database. The so- 
lutions typically have errors of < 0.03 mag in both passbands 
and negligible color terms. 

Table [2] summarizes the type of photometric solution for 
each mosaic. The table gives in Col. 1 the field name, in Col. 2 
the passband, in Col. 3 the number of default solutions, in gen- 
eral corresponding to non-photometric nights, and in Cols. 4-6 
the number of nights with photometric solutions derived from 
1- to 3- parameter fits and in Col. 7 the total number of nights 
used to observe pointings in a given mosaic. In most of the mo- 
saics solutions with at least 2 parameters are available thereby 
assuring an adequate calibration of at least some frames in 
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Table 2. Type of photometric solutions available for each mo- 
saic. 



Field 


Passband 


Default 


1-par 


2-par 


3-par 


total 


Deep la 


J 





6 





1 


7 


Deep lb 


.1 


2 


4 








6 


Deep2b 


.1 


1 








7 


8 


Deep2c 


.1 








4 


4 


8 


Deep3a 


J 





2 


4 


4 


10 


Deep3b 


J 





1 


4 


1 


6 


Shallow 


Deep la 


K s 





1 


2 





3 


Deep lb 




1 





3 


1 


5 


Deep2a 


^ 


2 


1 


2 


1 


6 


Deep2b 


K s 


1 


1 


4 


6 


12* 


Deep2c 


K s 


1 


1 


6 


6 


14* 


Deep3a 


K s 


2 


4 


5 


3 


14* 


Deep3b 


K s 


1 


4 


4 


1 


10* 


Deep 


Deep la 


K s 


1 


5 


2 


2 


10 


Deep lb 


Ks 


2 


1 


4 


3 


10 


Deep2b 


Ks 


1 


1 


4 


6 


12 


Deep2c 


Ks 


1 





3 


6 


10 


Deep3a 


Ks 


2 


3 


4 


3 


12 


Deep3b 


K s 


1 


3 


4 


1 


9 



includes Deep 



Table 3. Average photometric solutions 



Passband 


Zp 


k 


color 


J 


23.12±0.09 


0.08±0.04 


-0.02±0.05 


K s 


22.41±0.08 


0.08±0.06 


-0.01±0.24 



each mosaic. The only exception is Deeplb in /, where only 
1 -parameter solutions were possible, leading to a less accurate 
calibration of these data. In only one of the remaining cases 
(Deepla K s shallow) no 3-parameter solution is available. 

Table |3] summarizes the average photometric solutions de- 
termined for the entire period of observations. For each pass- 
band the table lists in Col. 1 the passband, in Col. 2 the mean 
and standard deviation of the measured zero points, in Col. 3 
the mean and standard deviation of the determined extinction 
coefficients, in Col. 4 the mean and standard deviation of the 
measured color terms. The listed standard deviations are a good 
indication of the uncertainty of the estimated zero points. 

The results of the science reduction are summarized in 
Table g] This table lists in Col. 1 the field name, in Col. 2 
the passband, in Col. 3 the number of pointings observed, in 
Col. 4 the total on-source time in ksec, in Col. 5 the number of 
nights in which the observations were carried out, in Col. 6 the 
number of raw exposures involved, and in Col. 7 the number of 
reduced images available. A total of 609 reduced images were 
produced from the 14567 science exposures corresponding to a 
total integration time of 874 ksecs (~ 243 hours). 

All the reduced images were inspected and graded and in 
Table [5] the results of the visual inspection are summarized 
grouping by mosaic. The table lists in Col. 1 the name of the 



Table 5. Quality assessment of reduced images. 



Field 


Passband 


A 


B 


C 


D 


Deepla 


J 


34 


6 








Deeplb 


J 


18 


12 





1 


Deep2b 


.1 


35 


6 








Deep2c 


.1 


12 


1 





1 


Deep3a 


J 


19 


2 


1 





Deep3b 


J 


28 


4 








Shallow 


Deepla 


K s 


16 


18 


4 


1 


Deeplb 


Ks 


10 


14 


6 


5 


Deep2a 


Ks 


28 


4 


5 





Deep2b 


Ks 


60 


11 


6 


1 * 


Deep2c 


Ks 


55 


6 


3 


2 * 


Deep3a 


Ks 


42 


10 


2 


* 


Deep3b 


K s 


38 


12 


5 


1 * 


Deep 


Deepla 


K s 


6 


14 


12 





Deeplb 


Ks 


17 


10 


3 


2 


Deep2b 


Ks 


41 


10 


6 


1 


Deep2c 


Ks 


26 


1 


1 





Deep3a 


Ks 


27 


7 








Deep3b 


K s 


25 


8 


3 


1 



* includes Deep 



field, in Col. 2 the passband, and in Cols. 3-6 the grades rang- 
ing from A (best) to D (worst). These grades reflect the overall 
cosmetic quality of the images including the background and 
the presence of other features that may affect the source extrac- 
tion. In addition to the grade a subjective comment is normally 
associated with the image, to describe any abnormality or fea- 
ture relevant to a user. Normally the comments are related to 
poor background estimation or remaining imprints of the fring- 
ing in particular in the K s -band images. Typically, the /-band 
images are of excellent quality with over 80% being graded 
A. The worse grades are usually due to the presence of bright 
vertical lines resulting from cross-talk effects associated with 
the presence of relatively bright stars in or near the field. In 
the Ar s -band ~ 65% are graded A and an additional ~ 20% are 
graded B. These grades, on average lower, reflect in particular 
that residual fringing is more common in this band. 

4. Final products 

4. 1. Images 

Final stacked images were produced from the nightly reduced 
images, described in the previous section, properly grouped in 
Stacking Blocks (SBs) as explained in detail in Paper I. The 178 
and 417 reduced images with grades better than D in the J- and 
ZTs-bands, were converted into 89 and 272 stacked (co-added) 
images. The process of producing the final stacks takes into ac- 
count differences in zero point that exist for individual images 
by correcting all frames to the zero airmass flux level. Non- 
photometric frames are scaled to the zero airmass flux level of 
the photometric frames, thereby assuring the photometric qual- 
ity of the produced images. 
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Field 


Passband 


# Pointings 


T (ksec) 


Nights 


Exposures 


# Images 


Deep la 


J 


16 


60.4 


7 


1006 


40 


Deep lb 


J 


13 


53.3 


6 


888 


31 


Deep2b 


.1 


16 


60.4 


8 


1007 


41 


Deep2c 


J 


12 


45.4 


8 


756 


14 


Deep3a 


J 


16 


63.4 


10 


1056 


22 


Deep3b 


J 


16 


57.5 


6 


958 


32 


Shallow 


Deep la 


K s 


36 


26.8 


3 


447 


39 


Deep lb 


A'v 


30 


20.7 


5 


345 


35 


Deep2a 


K s 


29 


20.0 


6 


334 


29 


Deep2b 


A'v 


36 


89.6 


12 


1494 


78* 


Deep2c 


K s 


49 


82.5 


14 


1375 


66* 


Deep3a 


K s 


36 


89.8 


14 


1498 


54* 


Deep3b 


K s 


35 


94.3 


10 


1572 


56* 


Deep 


Deep la 


K s 


14 


51.6 


10 


860 


32 


Deep lb 


A'. 


14 


52.5 


10 


876 


32 


Deep2b 


Av 


16 


75.7 


12 


1262 


58 


Deep2c 


A'v 


12 


57.3 


10 


956 


28 


Deep3a 


A, 


16 


75.9 


12 


1266 


34 


Deep3b 


K s 


16 


80.9 


9 


1349 


37 


* includes Deep 



Table 6. Quality assessment of stacked images. 



Field 


Passband 


A 


B 


C 


D 


Deep la 


J 


14 


2 








Deep lb 


.1 


11 


2 








Deep2b 


J 


16 











Deep2c 


.1 


11 


1 








Deep3a 


./ 


14 


2 








Deep3b 


/ 


14 


2 








Shallow 


Deep la 


A, 


17 


10 


8 





Deep lb 


As 


10 


11 


5 


1 


Deep2a 


As 


27 


1 


1 





Deep2b 


As 


33 


3 





* 


Deep2c 


As 


40 


3 


3 


* 


Deep3a 


As 


23 


8 


3 


2 * 


Deep3b 


K s 


26 


4 


3 


2 * 


Deep 


Deep la 


K s 


2 


5 


6 


1 


Deep lb 


Ks 


10 


4 








Deep2b 


Ks 


15 


1 








Deep2c 


Ks 


11 





1 





Deep3a 


Ks 


8 


6 


2 





Deep3b 


K s 


12 


2 


2 






includes Deep 



As before all images were visually inspected and graded. 
Out of the 361 stacked images covering the 7 DPS-NIR mo- 
saics, 268 were graded A, 56 B, 31 C, and 6 D. In Table[6]the 
grade distribution for each mosaic is given. The table lists in 
Col. 1 the field name, in Col. 2 the passband, and in Cols. 3-6 
the number of images in each mosaic with grade A-D, respec- 



tively. As in the case of the reduced images the comments made 
for the stacks are normally related to residual fringing and poor 
background subtraction, both of which are most frequent in the 
AVband. In /-band the most frequent problem is cross-talk fea- 
tures. 

The main attributes of the stacks produced for each mosaic 
are summarized in Table [7] The table gives in Col. 1 the field 
name, in Col. 2 the passband, in Col. 3 the number of pointings 
available with grade better than D, in Cols. 4 and 5 the aver- 
age and standard deviation of the FWHM in arcseconds and 
anisotropy (PSF rms) of the point-spread function measured in 
the final stacks, in Col. 6 the average and standard deviation 
of the limiting magnitude, m/, m , estimated for the final image 
stack for a 2" aperture, 5o detection limit in the Vega system, 
in Col. 7 the average and standard deviation of the complete- 
ness expressed as the fraction (in percentage) of observing time 
relative to that originally planned. 

The overall image quality is high with the average PSF 
having a FWHM of < 1" and an anisotropy typically less 
than 10%. Two of the shallow mosaics (Deepla and Deep2c) 
show a significantly larger FWHM. This is because in these 
mosaics there are a few frames taken in poor seeing condi- 
tions, and does not reflect the quality of these mosaics as a 
whole. The overall median limiting magnitudes in the Vega 
system are 22.16 mag (Jab = 23.06) in /-band; A^=19.57 mag 
(Ksab = 21 .41) for the shallow strategy; A^—20.32 mag (K s ab = 
22.16) for the deep strategy. This means that the proposed 
depth of K s ab =21.3 mag for the AVband shallow has on aver- 
age been reached, while the deep JK S surveys are about 0.3-0.5 
magnitudes shallower than originally planned (Jab = 23.4 and 
Ksab = 22.7 mag). 
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Table 7. Stacked image attributes and completeness for each mosaic only considering pointings with grade better then D. 



Field 


Passband 


$ Pointings 


PSF FWHM 


PSF rms 


" l lim \ vt S a / 


Completeness 


Deep la 


J 


16 


0.676±0.094 


0.086±0.025 


22.17± 0.23 


104%± 27% 


Deep lb 


.1 


13 


0.935±0.238 


0.073±0.034 


22.14± 0.56 


86%± 38% 


Deep2b 


.1 


16 


0.791±0.191 


0.110±0.041 


22.06± 0.16 


103%± 12% 


Deep2c 


.1 


12 


0.837±0.311 


0.090±0.034 


22.01±0.17 


96%± 11% 


Deep3a 


J 


16 


0.751±0.141 


0.090±0.033 


22.02± 0.49 


110%±40% 


Deep3b 


J 


16 


0.712±0.066 


0.076±0.016 


22.05± 0.17 


99%± 0% 


Shallow 


Deep la 


K s 


35 


1.275±0.466 


0.053±0.021 


19.57± 0.16 


103%± 25% 


Deep lb 


K s 


26 


0.890±0.198 


0.065±0.040 


19.38± 0.29 


97%± 30% 


Deep2a 


K s 


29 


0.749±0.274 


0.066±0.023 


19.43± 0.20 


95%± 11% 


Deep2b 


Ks 


36 


0.715±0.181 


0.064±0.025 


19.47± 0.16 


96%± 8%* 


Deep2c 


Ks 


46 


1.387±0.602 


0.057±0.023 


19.62± 0.14 


96%± 21%* 


Deep3a 


K x 


34 


0.737±0.092 


0.059±0.011 


19.39± 0.08 


100%± 0%* 


Deep3b 


Ks 


33 


0.688±0.097 


0.052±0.009 


19.30± 0.09 


100%± 0%* 


Deep 


Deep la 


K s 


13 


0.712±0.090 


0.095±0.026 


20.07± 0.23 


77%± 16% 


Deep lb 


K s 


14 


0.911 ±0.209 


0.089±0.026 


20.24± 0.31 


77%± 21% 


Deep2b 


K s 


16 


0.842±0.151 


0.144±0.039 


20.44± 0.17 


97%± 5% 


Deep2c 


K s 


12 


0.671 ±0.080 


0.133±0.043 


20.30± 0.20 


99%± 23% 


Deep3a 


Ks 


16 


0.712±0.159 


0.102±0.027 


20.34± 0.14 


98%± 9% 


Deep3b 


Ks 


16 


0.679±0.080 


0.104±0.020 


20.26± 0.12 


101%± 18% 


* includes Deep 



The completeness of the survey is measured not only in 
terms of depth but also in areal coverage. The latter can be as- 
sessed from Figs. IA.llA.3l in the appendix. These figures show 
the final coverage for each region and strategy, after discarding 
images with grade D. Note that in the case of Deepl the K s deep 
observations were done independently of those conducted for 
/fj-shallow. As can be seen the areal completeness and connec- 
tivity of the mosaics are high, with the exception of the shallow 
observations of Deepla and b, which are 70% and 53% areal 
complete, respectively. 



18.94 ±0.31 mag, with a typical coverage of ~ 25.8 square 
arcmin. 

Inspection of the CLASS_STAR parameter as a function of 
magnitude indicates that stars and galaxies can be adequately 
classified by using CLASS_STAR with a threshold of 0.90 at 
magnitudes K s < 17.0 in both the deep and shallow data sets. 
For the /-band we identify stars brighter than J = 18.0. At 
fainter magnitudes we take all objects to be galaxies. In the next 
section, these catalogs are used to assess the scientific quality 
of the survey data. 



4.2. Source lists 

Catalogs were extracted from the stacked images using the EIS 
Data Reduction System adopting the same procedures as de- 
scribed in Papers I and II. A description of the catalog format 
is found in the appendices of Paper I. These catalogs include 
only objects with S/N > 3 as computed from the magnitude er- 
ror. Table[8]summarizes some of the main characteristics of the 
extracted catalogs. The table gives in Col. 1 the field name, in 
Col. 2 the passband, in Col. 3 the average and standard devia- 
tion of the effective area in square arcmin, in Col. 4 the average 
and standard deviation of the number of objects, and in Col. 5 
the average and standard deviation of the 80% completeness 
limit in the Vega system. These completeness limits are for ex- 
tended objects (see Paper I for details). Their distribution for 
each strategy is illustrated in Fig. ^ There is a total of 89 / 
and 87 K s catalogs with mean values of the 80% completeness 
magnitude of 21.49±0.31 and 19.57±0.28 mag in the Vega 
system, with each catalog covering an area of ~ 28.8 square 
arcmin. For the shallow K s strategy there are 179 images with 
a mean 80% completeness magnitude in the Vega system of 



5. Discussion 

5. 1. Internal photometric comparison 

The consistency of the photometric calibration of the individual 
pointings, and therefore of the mosaic, was assessed by com- 
paring the magnitudes of objects in common between adjacent 
frames. For each mosaic the mean and standard deviation of 
the magnitude difference of all pairs of sources in common in 
adjacent frames were computed, considering only objects with 
a magnitude error < 0.2 (~ 5a). Ideally, one would use only 
stars for this, but because the number of objects in the over- 
laps is small, it proved necessary to include galaxies as well. 
In Table|5]these mean offsets and standard deviations are listed 
for each mosaic. The table lists in Col. 1 the field name, in 
Col. 2 the passband, in Col. 3 the number of overlaps in the 
mosaic, in Cols. 4 and 5 the mean and standard deviation of 
the magnitude offsets. From the table one finds that on average 
the offsets are small, typically < 0.05 mag. The scatter is typ- 
ically less than 0.1 mag. In all cases the mean offset is smaller 
than the standard deviation. Combined these results suggest a 
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Fig. 1. Distribution of the 80% completeness limiting Vega-magnitudes for each of the three observing strategies as indicated in 
each panel. 

Table 8. Summary of the properties of the catalogs extracted for each mosaic and band. 



Field 


Passband 


#Catalogs 


Eff. area 


#Objects 


80% compl. 


Deep la 


J 


16 


29.4± 1.0 


908± 110 


21.50± 0.23 


Deep lb 


.1 


13 


28.7± 1.8 


665± 145 


21.57± 0.51 


Deep2b 


.1 


16 


29.4± 1.0 


828± 90 


21.48± 0.23 


Deep2c 


.1 


12 


28. 8± 2.1 


81 1± 208 


21.41±0.25 


Deep3a 


.1 


16 


27.2± 2.2 


934± 178 


21.52± 0.38 


Deep3b 


J 


16 


29.4± 0.7 


907± 80 


21.44± 0.25 


Shallow 


Deep la 


K s 


35 


26.0± 1.1 


361± 161 


18.94± 0.38 


Deep lb 


K s 


26 


24.6± 1.7 


484± 112 


18.95± 0.24 


Deep2a 


K s 


29 


26.1± 1.6 


440± 91 


19.05± 0.22 


Deep2b 


K x 


20 


26.1± 1.2 


479± 78 


19.05±0.17 


Deep2c 


Ks 


34 


25.7± 1.6 


318± 188 


18.79± 0.47 


Deep3a 


K s 


18 


26.3± 0.5 


593± 98 


19.01±0.09 


Deep3b 


K s 


17 


26.7± 0.6 


551±53 


18.91±0.08 


Deep 


Deep la 


K s 


13 


29.4± 1.1 


996± 133 


19.40± 0.22 


Deep lb 


K s 


14 


28.9± 1.2 


660± 117 


19.57± 0.48 


Deep2b 


K s 


16 


28.1±0.8 


697± 61 


19.70± 0.21 


Deep2c 


K s 


12 


28.6± 1.7 


699± 80 


19.52± 0.21 


Deep3a 


Ks 


16 


26.9± 1.9 


804± 126 


19.63± 0.18 


Deep3b 


K s 


16 


29.8± 0.6 


932± 161 


19.63± 0.12 



* includes Deep 



homogeneous photometric calibration across the mosaics with 
fairly small field-to-field variations. 



In an attempt to improve the match among different zero 
points obtained on different nights, correct for data taken un- 
der non-photometric conditions, and he nce minimize field-to- 
field variations, the method employed by Madd oxet alJ (Q990) 
was utilized. Taking into account the computed corrections, one 
finds that, in general, the mean offsets do not change signif- 
icantly and the scatter decreases only marginally. This result 
may be partly due to the large uncertainty in each individual 
offset caused by the small number of sources in the overlap re- 
gions and the necessity to include galaxies. However, it may 
also indicate a good match among the original zero points. 



5.2. External comparison 

In order to externally verify the survey products, the source 
lists extracted from the final images as described in the previ- 
ous section were cross-correlated with the point source cata - 
log available from the 2MASS survey ( Skrutski eet all E006). 
This was carried out adopting a search radius of 1". A total of 
1056, 966, 1539 objects were found in common for the J, and 
the deep and shallow K s -band final images, respectively, over 
all regions covered by the survey. The differences in position 
and magnitude for objects in common and considered "good" 
in both catalogs (5a detections and in the case of 2MASS not 
contaminated by neighboring objects) were computed and their 
distribution as a function of right ascension and declination 
(left panels) and magnitude (right panels) are shown in Fig.|2] 
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Fig. 2. Astrometric (left column) and photometric (right column) comparison between objects in common between the EIS and 
2MASS surveys. The offsets are computed as 2MASS — EIS. The rows correspond to each adopted strategy as indicated. Dotted 
lines indicate the zero offset location, while the dashed lines in the right panels depict the mean offset computed for the bright 
objects as discussed in the text. In this plot magnitudes are in the Vega system. 



The mean offset in position is negligible and amounts to at most 
~ 0"07 for the shallow K s strategy. The scatter is of 0'.'23 and 
Of 27 in right ascension and declination, respectively, for all 
strategies. This is consistent with an astrometric accuracy of 
~ 0" 19 for the EIS images, limited by the typical internal ac- 
curacy of the reference catalog used. 



Considering objects brighter than / = 16.0 and K s = 14.5 
(Vega system) the measured magnitudes of the two data sets 
are in good agreement. At fainter magnitudes the scatter in- 
creases and the effects of the Malmquist bias introduced by 
the depth of the 2MASS data are seen. Objects with magni- 
tudes brighter than 12 mag are saturated in the EIS data and 
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Table 9. Mean magnitude offsets and standard deviations for 
each mosaic 



Field 


Passband 


# overlaps 


Mag. offset 


Std.dev. 


Deep la 


J 


24 


-0.02 


0.11 


Deep lb 


J 


15 


0.05 


0.11 


Deep2b 


.1 


24 


0.03 


0.05 


Deep2c 


.1 


12 


-0.03 


0.08 


Deep3a 


.1 


24 


-0.02 


0.06 


Deep3b 


J 


24 


0.02 


0.04 


Shallow 


Deep la 


K s 


36 


0.07 


0.13 


Deep lb 




21 


0.02 


0.18 


Deep2a 


^ 


39 


0.08 


0.12 


Deep2b 


K s 


8 


0.06 


0.22 


Deep2c 


K s 


34 


0.05 


0.13 


Deep3a 


Ks 


15 


0.06 


0.08 


Deep3b 


K s 


14 


0.01 


0.10 


Deep 


Deep la 


K s 


14 


0.00 


0.10 


Deep lb 


Ks 


16 


0.03 


0.11 


Deep2b 


Ks 


24 


0.03 


0.09 


Deep2c 


K s 


12 


0.00 


0.08 


Deep3a 


Ks 


19 


0.03 


0.07 


Deep3b 


K s 


24 


0.05 


0.08 



therefore discarded from the plots. For the deep strategies the 
systematic offset in magnitude is (discarding the outlying ob- 
jects with offsets in excess of 0.5 mag) ~ —0.04 mag with a 
scatter of 0.07 in K s and ~ 0.0 with a scatter of 0.08 in J, con- 
sistent with the estimated uncertainty of the overall zero-point 
(Sect. |3- For the shallow K s -survey one finds a mean offset 
of ^0.08 mag and a scatter of ~ 0.13 mag. The amplitude of 
the scatter is remarkably small, considering the values obtained 
in the previous section, and demonstrates the uniformity of the 
photometric calibration of the present data across the sky. 
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Fig. 4. Comparison of galaxy number counts in the J- (top) and 
Xj-band (bottom) between th e present work (solid c ircles with 
error bars); di amonds denote Vaisanenet al. (2000), triangles 
Martini (2001) and sauares llovino et alJ J2005I) . 



5.3. Comparison of galaxy and star counts 

An alternative way to evaluate the data is to compare the counts 
of galaxies and stars with those computed by other authors 
and/or predicted by models. For galaxies this is illustrated 
in Fig. The fi gure shows for each passband and strategy 
(columns) the galaxy counts obtained for each pointing (thin 
lines) grouped by mosaic (rows), as indicated in the panels 
of the leftmost panel of each r ow. These are com pared with 
the counts recently reported by llovino et al.l d2005l solid cir- 
cles). Here galaxies are 3o detections with CLASS_STAR< 
0.9. From the figure one can easily: 1) identify individual cases 
that depart significantly from the mean; 2) have a clear mea- 
sure of the variance introduced from field-to-field variations 
due to the relatively small field of view of SOFI; and 3) iden- 
tify the variation in depth of a few pointings, such as in the 
case of Deepla ^-shallow among others. Moreover, the rela- 
tively small scatter at intermediate magnitudes among the sub- 
fields forming a mosaic indicates that the relative photome- 
try between them is reasonable, reinforcing the evidence men- 
tioned above. More importantly, on averag e, there is a n excel- 
lent agreement with the counts of llovino et alJ d2005t) . which 



reach a comparable depth to the deep JK S strategy. This can be 
seen in Fig. |4] which shows the average counts obtained taking 
into account all fields and regions for the J- and /T v -band data. 
In all cases each field contributes only at magnitudes brighter 
than their respective 80% completeness limit, and no complete- 
ness corrections have been attempted. These mean counts show 
an excellent agreement with those recently obtained by differ- 
ent au thors JVaisanen et all 120001 iMartiniL 1200 it llovino et all 
2005), being further strong evidence of the overall quality of 
the data. The mean counts for the J- and A^-bands are given 
in Tables [Tolandll lllisting in Col. 1 the center of the magni- 
tude bin in the Vega system, in Col. 2 the number of galaxies 
per 0.5 mag and per square degree, and in Col. 3 the standard 
deviation among the contributing fields. 

Finally, the corresponding stellar counts, sources with 
CLASS_STAR> 0.9, are shown in Fig. \5\ The counts (solid 
circles connected by a solid line) are compared to pre dic- 
tions adopting the model proposed by Gira rdi et all (|2005, thin 
lines). Despite the large scatter due to the small field-of-view, 
the agreement is remarkable, considering the different regions 
and that the model was parameterized using optical data. 
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Fig. 3. Galaxy counts for each individual pointing grouped by mosaic (lines) compared with those of llovino et alJ J2005l solid 
circles). 



6. Summary 

This paper presents the infrared data accumulated by the 
Deep Public Survey conducted by the EIS program. The pa- 
per presents the observations, reductions, and some of the 
main products generated and administrated by the EIS Data 
Reduction System which allows the unsupervised reduction 



of large amounts of data. For instance, the data set described 
in the present paper was reduced overnight using eight dual- 
processors computers running Linux. The highthroughput of 
the system also makes it ideal for handling the new generation 
of infrared surveys using wide-field cameras. In general, the 
depth of the survey is relatively homogeneous and the cover- 
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Fig. 5. Comparison of star counts (filled circles connected by solid line) for each mosaic (rows) and strategy (columns) compared 
with model predictions bv lGirardietaD ( l2005t thin line). 



age rather uniform with exception of the shallow coverage of 
Deepl. However, considerable improvement in coverage and 
photometric accuracy could have been reached if the survey 
was carried out in service rather than visitor mode allowing a 
better optimization depending on the actual weather conditions, 



in particular ensuring sufficient photometric data for reliable 
calibration of all the fields. 

The near-infrared dataset being released is one of largest 
currently available and its quality has been assessed qualita- 
tively by visual inspection of the individual images and quanti- 
tatively by direct comparison with measurements made by the 
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Table 10. The derived raw galaxy counts for the /-band. The 
counts are given per mag per square degree. Each field only 
contributes to the bins brighter than their 80% completeness 
limit. No completeness correction has been attempted. 





NT 


Gpj 


in 0^ 
1U.Z0 


u.u 


u.u 


1 n 7^ 


u.u 


U.U 


1 1 

1 1 .Zj 


u.u 


U.U 


1 1 7^ 


u.u 


U.U 


1 Z.ZJ 


Q 

z.y 


ZO.7 


1 1 7^ 


Z.ft 


zo.u 


1 J.Z.J 


Z.ft 


zo.u 


1 ^ 7S 


1 Q 


jU.o 


1 A T\ 


I 1 1 

I I .Z 


JZ.Z 


1 A 7^ 
14. / 


47 A 


1 in Q 
1 lU.o 


1 J.ZJ 


13.1 


Ijj.j 


1 ^ 7^ 

ID. ID 


7 


100. 4 


10. ZO 


Zoj.j 


Zj J.O 


1 ^ 7S 

10. / J 


.J 


jOI.U 


17.25 


766.0 


463.1 


17.75 


1350.6 


606.4 


18.25 


2059.3 


746.3 


18.75 


3189.6 


1083.7 


19.25 


4988.9 


1731.7 


19.75 


7804.7 


2398.9 


20.25 


11671.5 


3354.3 


20.75 


16064.9 


3381.6 


21.25 


18383.9 


5991.3 


21.75 


12933.6 


8897.6 


22.25 


10949.9 


10678.0 



Table 11. The derived raw galaxy counts for the /T s -band. The 
counts are given per mag per square degree. Each field only 
contributes to the bins brighter than their 80% completeness 
limit. No completeness correction has been attempted. 



Mag 


N 




11.25 


0.9 


14.6 


11.75 


4.1 


33.6 


12.25 


2.1 


24.2 


12.75 


11.9 


54.9 


13.25 


20.2 


77.1 


13.75 


51.8 


118.7 


14.25 


121.2 


209.2 


14.75 


224.4 


254.5 


15.25 


471.6 


421.6 


15.75 


747.9 


509.5 


16.25 


1282.1 


664.6 


16.75 


2144.8 


932.2 


17.25 


3598.5 


1366.0 


17.75 


5721.6 


1832.1 


18.25 


8872.0 


2904.3 


18.75 


11606.3 


3772.5 


19.25 


10333.8 


5998.6 


19.75 


8330.2 


5731.2 


20.25 


5295.8 


5289.0 



2MASS survey limited to the bright end and by the comparison 
of galaxy and star counts with those of other authors and with 
model predictions. 

The depth and areal coverage in the near-infrared com- 
bined with the deep multi-wavelength optical data reported in 
Paper III make this data set extremely u seful to further exp lore 
the nature of extremely red objects (e.g. iKong et all l2006t) and 
to search for distant galaxy clusters. With this paper, all the in- 
frared data accumulated by the EIS project up to late 2004 are 
now in the public domain 2 . 
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Appendix A: Coverage per region 

In this appendix, figures giving an overview of the final coverage for 
each mosaic are presented. Each figure covers one region (Deepl, 
Deep2 or Deep3). For Deepl there are three figures (J, deep K s , and 
shallow K s \ while for Deep2 and Deep3 only two figures (J and K s ) 
are available since here the deep K s strategy fills the central parts of 
the shallow mosaics. The relative sizes among the fields reflect their 
relative sky-coverage. 

From the figures it is clear that the areal completeness of the mo- 
saics are in general high, even though the Deepl shallow K s mosaics 
suffer from incompleteness, leading to large difficulties in creating 
mosaics and carrying out an internal re-calibration to establish a com- 
mon magnitude system for all frames. The centers of each subfield can 
be found in Table Q 
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Fig. A.2. The coverage in the Deep2 region (north up, east left). The J coverage (52'.5 x 19'3) is on the top and K s (90'8 x 3 1 '.2) 
on the bottom. The fields in the upper panel are Deep2b (left) and Deep2c (right). The fields in the lower panel are Deep2a, 
Deep2b and Deep2c from left to right. In the ^-band the central parts in K s for the mosaics Deep2b and Deep2c are the deep 
strategy, for Deep2a only shallow data exists. 
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Fig.A.3. The coverage in the Deep3 region (north up, east left) with J (52'.0 x 19'.1) on the top and K s (60 'A x 27 '3) on the 
bottom. The fields are Deep3a (left) and Deep3b (right). For the 7T. s -band the central parts correspond to the deep strategy. 



